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Numerical evidence for Goldstone-mode effects in the three-dimensional 
0(4)-model* 

Jurgen Engels and Tereza Mendes a 

a Fakultat fur Physik, Universitat Bielefeld, D-33615 Bielefeld, Germany 

We investigate the three-dimensional 0(4)-model on 24 3 -96 3 lattices as a function of the magnetic field H. 
Below the critical point, in the broken phase, we confirm explicitly the H 1 ^ 2 dependence of the magnetization and 
the corresponding H' 1 ^ 2 divergence of the longitudinal susceptibility, which are predicted from the existence of 
Goldstone modes. At the critical point the magnetization follows the expected finite-size-scaling behavior, with 
critical exponent 5 = 4.87(1). 



The continuous symmetry present in the O(N) 
spin models 
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— where Si are unit vectors in an iV-dimensional 
sphere — gives rise to the so-called spin waves: 
slowly varying (long-wavelength) spin configura- 
tions, whose energy may be arbitrarily close to 
the ground-state energy. In two dimensions these 
modes are responsible for the absence of sponta- 
neous magnetization, whereas in d > 2 they are 
the Goldstone modes associated with the spon- 
taneous breaking of the rotational symmetry for 
temperatures below the critical value T c . Due to 
the presence of Goldstone modes, a diverging sus- 
ceptibility is induced in the limit of small external 
magnetic field H for all T < T c , i.e. everywhere on 
the coexistence curve (see for example Q). More 
precisely: not only the transverse susceptibil- 
ity, which is directly related to the fluctuation of 
the transverse (Goldstone) modes, diverges when 
H ->0 



X t(T<T c ,H) 



M (T, H) 
H 



= 0{H- r ) 



but also the longitudinal susceptibility, given 
in terms of the magnetization M by \l = 
dM/dH , is predicted to diverge for 2 < d < 4. 
In d = 3 the predicted divergence is H^ 1 / 2 , or 
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cquivalently, the behavior for the magnetization 
will include a "singular" term of order H 1 / 2 

M(T < T C ,H) = M(T,0) + cH 1/2 . 

Indication of this behavior was seen in early nu- 
merical simulations of the three-dimensional 0(3) 
model at low temperatures M . 

Here we verify explicitly the predicted singu- 
larities — or Goldstone-mode effects — for the 
three-dimensional 0(4) case, by simulating the 
model in the presence of an external magnetic 
field and close to the critical temperature T c . The 
low-temperature singularities are still present at 
T <T C . We are able to compute magnetic crit- 
ical exponents directly, by varying H at T c , and 
we also use the observed Goldstone-effect behav- 
ior to extrapolate our data to H — > and ob- 
tain the zero-field critical exponents in very good 
agreement with the existing values. Using these 
exponents, we verify finite-size scaling at T c . Our 
motivation for considering the three-dimensional 
0(4) model is that it is expected to be in the same 
universality class as QCD at finite temperature 
for two degenerate quark flavors J3j , the magnetic 
field in the spin model corresponding to the quark 
mass in QCD. 

Our simulations are done using the cluster al- 
gorithm. We compute the magnetization along 
the direction of the magnetic field as 

M = y <^S i -H> . 
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Figure 1. Magnetization for T <T C . 

We note that, due to the presence of a nonzero 
field, the magnetization defined above is nonzero 
on finite lattices, contrary to what happens for 
simulations at iJ = 0, where one is led to consider 
< I 2 I >i which overestimates the true 
magnetization. 

The susceptibilities are obtained from the spin 
correlation functions 

g l {x) = <s\s\> - <sl><sl> 

Gt{x) = 1<S^-S^> 

by xl = dM/dH = Gl(0); xt = G T (0) . We 
use for T c the value 

1/T C = J c = 0.93590, 

obtained in simulations of the zero-field model ||] . 

In Fig. [l] we show our data for the magnetiza- 
tion for temperatures T <T C plotted versus H 1 / 2 . 
Inverse temperatures are given by J = J c , 0.95, 
0.98, 1.0, 1.2, starting from the lowest curve. We 
have simulated at increasingly larger values of L 
at fixed values of J and H in order to eliminate 
finite-size effects. Except for the curve at J = 1.2 
only these largest-lattice values are shown. Sym- 
bols are as explained in Fig. |[ 
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Figure 2. Finite-size-scaling plot for the magne- 
tization at T c . 

The finite-size effects for small H do not dis- 
appear as one moves away from T c , but rather 
increase. Solid lines connect the y-axis to the last 
point (a filled square) included in our fits of the 
data for small H . It can clearly be seen that the 
predicted behavior (linear in H 1 / 2 ) holds close 
to H = for all temperatures T < T c consid- 
ered. We thus see evidence that the Goldstone- 
mode effects are observable even very close to T c . 
The straight-line segments become shorter as T c 
is approached from below, and at T c the magne- 
tization vanishes as a power of H, as expected. 

We have fitted the data from the largest lattice 
sizes at T c to the scaling behavior 

M(T C ,H) ~ H 1 / 6 , 

obtaining the exponent S — 4.87(1), in agreement 
with Q. The straight-line fits shown in Fig. |l| 
are used to extrapolate our data to the zero-field 
limit, and the corresponding zero-field scaling law 

M(T<T C ,H = 0) - (T c -Tf 

yields j3 — 0.38(1), also in agreement with Q. 
The finite-size effects at T c are very well described 
by the finite-size-scaling ansatz 

M(T C ,H;L) = L-( i ' v Q M {HL^I v ) , 
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Figure 3. Susceptibilities \t and Xl at low T. 

as shown in Fig. |2|. We see no corrections to finite- 
size scaling. In Fig. ^ we show a typical behavior 
for xt and xl at low temperature, for J — 0.98. 
Symbols are as in Fig. |. The lines M/H (a factor 
3 is included for clarity) and dM / dH (from the fit 
in Fig. are shown, for comparison respectively 
with xt and xl ■ Both lines are obtained from our 
"infinite-volume" data. Notice the large finite- 
size effects in xl- 

Finally, we verify numerically the claim in 
Ref. |U, that the Goldstone-mode singularity at 
low temperatures is consistent with the magnetic 
equation of state 

where y = H/M s . In Figs. [| and || we show 
respectively the behavior of xl and the combina- 
tion above for J = 0.95. We see that the two per- 
turbative predictions are well satisfied at T < T c . 
Note that the line in Fig. || is not obtained using 
the perturbative coefficients given in fjj , but from 
a fit of our data. This idea of finding nonpertur- 
bative coefficients from fits to perturbative forms 
is very useful in obtaining an expression for the 
scaling function for this model ||. 

This work was supported by the DFG under 
Grant No. Ka 1198/4-1 and by the TMR-Network 
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Figure 4. The susceptibility xl at J = 0.95. 
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Figure 5. (I/^M'-^l at J = 0.95. 
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